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Coherent Synchrotron Radiation (CSR) can play an important role by not only increasing the 
energy spread and emittance of a beam, but also leading to a potential instability. Previous studies 
of the CSR induced longitudinal instability were carried out for the CSR impedance due to dipole 
magnets. However, many storage rings include long wigglers where a large fraction of the synchrotron 
radiation is emitted. This includes high-luminosity factories such as DAPHNE, PEP-II, KEK-B, 
and CESR-C as well as the damping rings of future linear colliders. In this paper, the instability 
due to the CSR impedance from a wiggler is studied assuming a large wiggler parameter K. The 
primary consideration is a low frequency microwave-like instability, which arises near the pipe cut-off 
frequency. Detailed results are presented on the growth rate and threshold for the damping rings of 
several linear collider designs. Finally, the optimization of the relative fraction of damping due to 
the wiggler systems is discussed for the damping rings. 
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I. INTRODUCTION 

When an electron bunch passes the curved trajectory 
inside a dipole, the coherent synchrotron radiation (CSR) 
induces an impedance 0,0,0)0- The CSR impedance 
is known to impact single pass bunch compressors where 
the beam currents are extremely high [a, IS HI but 
it is also possible that CSR might cause a microwave- 
like beam instability in storage rings. A theory of such 
an instability in a storage ring has been recently pro- 
posed in Ref. |9| where the impedance is generated by 
the synchrotron radiation of the beam in the storage ring 
bending magnets. 

A similar instability may arise due to coherent syn- 
chrotron radiation in long wigglers. Many storage rings 
have used large wiggler systems to reduce the damping 
times and control the beam emittance. In particular, the 
high luminosity colliding beam factories DAPHNE [Icj . 
PEP-II H2, KEK-B Gj, CESR-C [13 all have wiggler 
systems that can radiate significantly more power than is 
radiated in the arc bending magnets. Similarly, damping 
ring designs for future linear colliders rely on large wig- 
gler systems to obtain the required damping rates where, 
again, more power can be radiated in the wiggler systems 
than in the arc bending magnets. 

In a previous study, the wakefield and impedance in a 
wiggler with large parameter K has been obtained [T^ 
based on the work of Saldin et al. 0. In this paper, 
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we study the impact of the wiggler synchrotron radiation 
impedance on the beam longitudinal dynamics in rings 
with dipoles and wigglers. We focus our attention on the 
limit of a large wiggler parameter K because this is the 
most interesting case for many applications. 

The paper is organized as follows. In Sec. [H] we will 
briefly review the theory of the beam instability in a ring 
due to the CSR impedance 0. Using the synchrotron ra- 
diation impedance of a wiggler we study the beam 
instability problem for a ring including both dipoles and 
wigglers in Sec. IIHI In Sec. II VI we then focus on a low- 
frequency microwave-like CSR instability near the pipe 
cut-off frequency which is the primary limitation. Fi- 
nally, in Sec. El we discuss effect of the CSR instability 
on the design optimization of damping rings for future 
linear colliders. Finally, we give some discussions and 
conclusion in Sec. IVII 



II. LONGITUDINAL BEAM INSTABILITY 

We consider the longitudinal dynamics of a thin coast- 
ing beam. The beam can be described with a longitudinal 
distribution function p(y, s, z). The positive direction for 
the internal coordinate s is pointing to the direction of 
the motion. The relative energy offset of a particle hav- 
ing energy E with respect to the nominal energy Eq is 
expressed as v = (E — Eq)/Eq. The position of the ref- 
erence particle in the beam line is z = c t with c equal to 
the speed of light in vacuum. 

The longitudinal beam instability due to the syn- 
chrotron radiation impedance in a ring has been studied 
theoretically by Stupakov and Heifets In this section, 
we briefly review and reproduce the equations we need 
for our study in this paper. The reader is advised to refer 
to reference for the details of the theory. 
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For a matched beam, one can use a 1-D Vlasov equa- 
tion to describe the evolution of the longitudinal distri- 
bution function p(y, s, z). 
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where, the slippage factor is defined as 
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with a equal to the momentum compaction factor. For 
ultrarelativistic beams, we have r/ ss a. Based on our def- 
inition of s, the sign convention for a is then the follow- 
ing. If a particle with a higher energy than the nominal 
energy, i.e., v > 0, would go to the head of the bunch, 
then the beam line provides a negative a. According to 
this convention, then in a simple bending magnet, we will 
have a > 0, while in a wiggler, we have a < 0. 

The distribution function p is written as a sum of 
the equilibrium distribution function and a perturbation, 
p = po(v) + pi{v, s, z), with pi <C po- We look for a per- 
turbation as 

^ — iuj zj c-\- ik s 



Pi = Pie 
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where k is the wavenumber, uj the frequency. Linearizing 
the 1-D Vlasov Eq. Q leads to the dispersion relation 



i tq c Z(k) 



dv- 



dpo I dv 
) + ckrjv ' 



(4) 



which determines the existence of a solution as in Eq. 
©. Here, r i=s 2.82 x 10~ 15 m is the classical electron 
radius and 7 is the Lorentz factor. In addition, 
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is the CSR impedance and the wake Green function w(s) 
describes the interaction of two particles due to the syn- 
chrotron radiation where w(s) 7^ for s > while 
w(s) — for s < 0. The positive values of w(s) cor- 
respond to the energy loss and the negative values imply 
the energy gain. 

We assume that the initial energy distribution function 
is a Gaussian, i.e., po = nct/(V2nvo) x exp(— v 2 /2 Vq), 
where no is the linear density, i.e., the number of particles 
per unit length. In this case, Eq. Q can be rewritten 
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where, A = 710 ^/(M 7 v^); Q, = ujj{ck \r]\ v§) and p — 
v/v n . The upper (lower) sign in Eq. 10 refers to the case 
of a positive (negative) r\. We refer the reader to Ref. 
for further details. 
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FIG. 1: The imaginary part of the normalized frequency fi as 
a function of the normalized wavenumber k/ko for the NLC 
main damping ring where ko is the on-axis wiggler funda- 
mental radiation wavenumber defined in Eq. 1141 . The solid 
curve includes the entire CSR impedance while the dotted 
and dashed curves include either the steady state dipole CSR 
impedance or the wiggler CSR impedance, respectively. The 
inset shows a blow up of the low frequency region where the 
beam is unstable. 



III. STORAGE RINGS WITH WIGGLERS 

As mentioned, a number of storage rings utilize long 
damping wigglers where a significant fraction of the en- 
ergy loss per turn is emitted. As concrete examples, 
we study the NLC main damping ring 0], the TESLA 
dam ping ring ^3] > and the KEK ATF prototype damping 
ring [13. Parameters are given in Table UJ 

In our paper, we use the steady state CSR impedance 
and assume a distributed impedance model. For a dipole, 
the steady state CSR impedance is Q, Q 



Z D (k) = -iA 



fcl/3 
R 2/3 



with 



A = 3-V*r (fj(V3i-l) 



The wiggler impedance Zw{k) is computed in Ref. [I 
Hence, the total impedance is then 

Z(k) = Z D (k)^ + Z w (k) L]l 



c 



(7) 

(8) 
-4j. 

(9) 



where, R, 0, L\y and C are the dipole bending radius, 
the total bending angle, the wiggler total length and the 
damping ring circumference given in Table [I] respectively. 

To study the instability, we numerically solve the dis- 
persion relation Eq. Since we are assuming a coast- 
ing beam model, we only consider instability wavelengths 
short compared to the bunch length. 
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NLC 


TESLA 


ATF 


Circumference C/km 


0.3 


17 


0.14 


Dipole radius R/m 


5.5 


80 


5.7 


Total bending angle 0/2-7T 


1 


5/3 


1 


Momentum compaction a/10 -4 


2.95 


1.2 


19 


Synchrotron frequency Q s /kHz 


3.5 


0.8 


17.4 


Extracted X emittance 7e :E /10 -6 m 


3 


8 


5 


Extracted Y emittance ■ye y /W~ s m 


2 


2 


5 


Energy E/Gev 


1.98 


5 


1.3 


Energy rms spread 4 


9.09 


9 


6 


Bunch rms length a z /mm 


3.6 


6 


5 


Particles in a bunch N e /10 10 


0.75 


2 


1 


Wiggler peak field B w /T 


2.15 


1.5 


1.88 


Wiggler period X w /m 


0.27 


0.4 


0.4 


Wiggler total length L w /m 


46.24 


432 


21.2 


tit- i n c i • n 1 

Wiggler p-function fj XlW /m 


1.87 


6.67 


6 


Pipe radius b/cm 


1.6 


2 


1.2 


Fw 


2.2 


13.4 


1.8 


Cut-off wavelength A c /mm 


4.9 


1.8 


3.1 


Threshold at cutoff (wiggler off) N t /10 10 


0.60 


27.44 


0.95 


Threshold at cutoff (wiggler on) N t /W 10 


0.52 


24.56 


0.76 


Growth time at cutoff (wigger off) r/fj.s 


54.9 


N/A 


34.3 


Growth time at cutoff (wigger on) r/fj,s 


32.9 


N/A 


6.5 



TABLE I: Parameters and results for the NLC main damping ring P5, the TESLA damping ring [12], and the KEK ATF 
prototype damping ring [l8|. The parameter F w is the ratio of the radiation power emitted in the wiggler to that emitted in 
the arc bending magnets defined in Eq. il llil . 



In Fig. the imaginary part of SI is plotted as a 
function of the instability wavenumber for the NLC main 
damping ring using the parameters listed in Table [I] At 
low frequencies, the dipole CSR impedance dominates 
while at shorter wavelengths the wiggler CSR impedance 
is usually more important. In the region where Im(f2)< 
0, the beam is stable and this is true for all regions except 
at the longest wavelengths as shown in the inset. This 
low frequency instability will be discussed in the following 
section. 

Similar calculations were made for the TESLA "dog- 
bone" damping ring, where the total bending angle is 
about 8 = 107r/3, and the KEK ATF prototype damping 
ring. For the design current in the TESLA damping ring, 
the impedance from the dipoles and wigglers will not 
drive an instability while the results for the ATF are very 
similar to those of the NLC damping ring shown in Fig. 
□ 

One interesting effect can be seen at the wiggler radia- 
tion fundamental frequency and the odd harmonics where 
the ring is actually stabilized by the wiggler impedance. 
This is opposite to the single-pass behaviour through a 
wiggler where there is an instability at the wiggler fun- 
damental frequency usually refered to as the FEL insta- 
bility. The effect will be discussed further in Appendix lAl 



but is a direct result of the positive sign of the momentum 
compaction in the ring. 

To illustrate this, Fig. [21 is a plot for the NLC damp- 
ing ring with identical wiggler and arc parameters but the 
momentum compaction is assumed to be opposite in sign. 
Here, one can see that at low frequencies the growth is 
slightly lower. In contrast, the wiggler impedance makes 
the system less stable at the wiggler fundamental unlike 
the case with positive momentum compaction. Further- 
more, if the magnitude of the momentum compaction is 
reduced, the system will become unstable similar to that 
illustrated in Appendix iBl This is similar to the FEL in- 
stability however, it is also noted in Appendix^that our 
theory does not correctly treat the full FEL instability. 

IV. LOW-FREQUENCY CSR INSTABILITY 

As seen in Figs. ^] and [21 the instability is most impor- 
tant at relatively low frequency. The longest instability 
wavelength that is possible is determined by two effects: 
first, the instability wavelength has to be short compared 
to the bunch length and, second, the vacuum chamber 
causes an exponential suppression of the synchrotron ra- 
diation at wavelengths A greater than the 'shielding cut- 
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FIG. 2: The imaginary part of the normalized frequency 
Q, as a function of the normalized wavenumber k/ko for 
the NLC main damping ring with negative momentum com- 
paction where ko is the on-axis wiggler fundamental radiation 
wavenumber defined in Eq. (I14H . The solid curve includes the 
entire CSR impedance while the dotted and dashed curves 
include either the steady state dipole CSR impedance or the 
wiggler CSR impedance, respectively. The inset shows a blow 
up of the low frequency region where the beam is unstable. 

off' 

A c < AV2b{b/R) 1/2 . (10) 

Here, R is the dipole bending radius, and b is the vac- 
uum chamber half height. The numerical coefficient of 
Eq. HlOfl assumes that the vacuum chamber is made up 
of two infinitely wide plates. Different cross sections give 
different numerical factors |23j. Given the previous dis- 
cussion, the threshold will be the lowest at the smaller of 
the bunch length or the 'shielding cutoff' wavelength. 

For the NLC main damping ring, we find that pertur- 
bations with wavelengths A > 3.5 mm are not stable due 
to the dipole CSR impedance alone. Adding the CSR 
impedance from the wiggler causes perturbations with 
wavelengths A > 2.6 mm to be unstable. In Fig. |3 the 
threshold particle number is plotted as a function of the 
perturbation wavelength. Next, in Fig. 0] the growth 
rate, which is denned as the inverse of the time needed 
for the perturbation to grow by a factor of e 

- = Im(n)cfc|i7|i/ . (11) 

T 

is plotted versus the perturbation wavelength. Based on 
the parameters in Table [I] and Eq. I|1U|) , the 'shielding 
cutoff' wavelength is computed to be A c ~ 4.9 mm. At 
this cutoff wavelength, the threshold currents and growth 
time are summarized in Table [I] The growth time is sig- 
nificantly faster than the synchrotron period, in agree- 
ment with the analysis for a microwave- like instability. 



FIG. 3: The threshold particle number as a function of the 
CSR wavelength for the NLC main damping ring. The dashed 
curve is the result for the dipoles only, while the solid curve 
takes into account of the contributions from the dipoles and 
the wigglers. The vertical straight line is the approximate cut- 
off wavelength according to Eq. 11111 . The horizontal straight 
line is the nominal number of particles per bunch: 0.75 x 10 10 . 

It is clearly seen that the threshold current decreases as 
we approach the longer wavelength perturbations, how- 
ever, the growth rate is not monotonic. This is the result 
of two opposite mechanisms: one is the energy modula- 
tion growth due to the CSR impedance, and the other is 
the Landau damping. For very long wavelength pertur- 
bation, Landau damping effect is small, hence, we could 
expand the denominator of the integrand in Eq. to 
get 

lo = ^ic 2 n b r r)kZ(k)/"f, (12) 

which indicates that r _1 cx Im(u;) cx \J kZ(k), i.e., the 
growth rate will increase for a shorter wavelength per- 
turbation, since Z(k) cx k e , with e > 0. This is shown 
in Fig. 0] when we approach from the long wavelength 
perturbation to the short wavelength perturbation, the 
growth rate increases. 

However, this process stops when Landau damping be- 
comes effective. The finite energy spread in the beam will 
produce a phase mixing due to the slippage and this will 
destroy density modulation due to the CSR induced en- 
ergy modulation. The Landau damping due to the phase 
mixing is more serious for short wavelength perturba- 
tions. This damping can be seen in the second term of 
Eq. Q , or more clearly, in the denominator of the dis- 
persion relation in Eq. J3J • This is demonstrated in Fig. 
0] the growth rate finally decreases when we approach 
very short wavelength, and eventually, the system be- 
comes stable. 

For the KEK ATF prototype damping ring, the cut-off 
wavelength would be about A c f=a 3.1 mm according to Eq. 
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FIG. 4: The growth rate as a function of the CSR wavelength 
for the NLC main damping ring. The dashed curve is the 
result for the dipoles only, while the solid curve takes into 
account of the contributions from the dipoles and the wigglers. 
The vertical straight line is the approximate cutoff wavelength 
according to Eq. 1101 . 



(CHI. In Fig. the threshold particle number is plotted 
as a function of the perturbation wavelength. Taking the 
dipole CSR impedance alone, for the single bunch charge, 
the instability sets in for perturbations with wavelengths 
A > 2.8 mm. Adding the wiggler CSR impedance, the 
electron beam would be unstable for perturbations with 
wavelengths A > 1.9 mm. Other results are summarized 
in Table IJ 

It is interesting to note that in both the NLC and the 
ATF damping rings roughly twice as much synchrotron 
radiation power is emitted in the wiggler as in the arc 
dipoles. However, the instability threshold is not dramat- 
ically impacted by the additional radiation in the wiggler 
and decreases by less than a factor of two in each case. 
This arises because of the very different low frequency 
behavior of the impedances. 

We can use a simple scaling analysis to understand 
this. According to Eq. JJJ, the dipole CSR impedance 
scales as Zu(k) oc fc 1 / 3 . Similarly, the low frequency 
behavior of the wiggler impedance is |l4j 



Zw(k) = 7r k v 



fen 



2i , ( k 

7T \k 



(13) 



which is accurate enough for practical calculation in k g 
[0, 0.1 fco], where ko is the wiggler fundamental radiation 
wavenumber: 



k a = 2j 2 k w /(l + K 2 /2) 



(14) 



and the wiggler parameter K is approximately K rs 
93.4 i?^ X w , with B w the peak magnetic field of the wig- 
gler in units of Tesla and A^, the period in meters. Thus, 



FIG. 5: The threshold particle number as a function of the 
CSR wavelength for the KEK ATF prototype damping ring. 
The dashed curve is the result for the dipoles only, while the 
solid curve takes into account of the contributions from the 
dipoles and the wigglers. The vertical straight line is the 
approximate cutoff wavelength according to Eq. 1101 . The 
horizontal straight line is the nominal number particles per 
bunch: 1.0 x 10 10 . 



the wiggler CSR impedance scales as Ke(Zw(k))<x k, and 
Im(Zn/(fc))oc fclog(fc), which have a weaker scaling than 
Zu(k) oc /c 1 / 3 when k — > 0, i.e., the CSR impedance from 
the wiggler is weaker than that from the dipole when we 
approach low frequency region. 

Hence, it suggests an optimization of the damping ring 
design where a larger fraction of radiation is emitted in 
the wigglers as will be discussed subsequently. 



V. DAMPING RING OPTIMIZATION 

Damping rings are used in linear colliders to attain 
the very small transverse emittances that required to at- 
tain the desired small spot sizes at the collision point. 
In the damping rings, the emission of synchrotron radia- 
tion damps the phase space volume of the injected beams 
towards the equilibrium values. An e- folding reduction 
of the transverse phase space occurs after the emission 
of synchrotron radiation energy comparable to the beam 
energy. 

There are usually three conflicting requirements on the 
damping rings: first, the equilibrium omittances must 
be small which usually requires weak bends and strong 
focusing; second, the dynamic aperture must be large to 
accept the injected beams and this is difficult to attain in 
a strong focusing storage ring; and, third, the damping 
rates must be fast so that the incoming phase space and 
transients can be damped before the beam are extracted 
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(15) 



and accelerated to the collision point. 

To achieve these conflicting requirements most damp- 
ing ring designs have been forced to include long damping 
wigglers where a significant fraction of the energy loss per 
turn is emitted. Of course, the ring characteristics can 
be greatly modified by changing the ratio of the damp- 
ing due to the wiggler versus that due to the bending 
magnets. However, historically the wiggler systems have 
been a concern because they are a significant source of 
non-linearity and because they generate intense radiation 
which may directly or indirectly impact the beam. 

Other work has improved the confidence in the dy- 
namic aperture predictions and has shown that properly 
designed wigglers should not impose a limitation |24j . 
Thus, we can further optimize the ring designs based on 
quantitative evaluations of the effect of the wiggler radi- 
ation fields. 

To study the ring optimization, we use expressions de- 
rived in Ref. [25j for the damping ring parameters. The 
primary constraint is the damping times. Because all 
linear collider designs presently only consider flat beams 
with a x 3> (Tj,, the tightest requirement is on the vertical 
damping time which, using Eq. (11) of Ref. p5j . can be 
written 

(2.89 x 10 12 kG)C |8j 
Ty ~ \B a \^{l + F w )c 2^ ' 

where B a is the arc bending field, C is the ring circum- 
ference, and 8 is the total bending angle of the arcs. 
In addition, F w is a parameter equal to the ratio of the 
damping due to the wiggler over the damping due to the 
arc bending magnets which can be written: 

F w = h w /h a ■ (16) 

Here, and l2 W are the second synchrotron integrals 
calculated over the arcs and the wigglers respectively, 
which can be written 

L B 2 

hw * 2(i# • (17) 

where {Bp) is the standard energy dependent magnetic 
rigidity, {Bp) ss 0.0177 (kG m), and B a , B Wl and L w are 
the magnetic field of the bending dipole, the peak wiggler 
field, and the wiggler length. 

According to Eq. Q15[l. we can scale the bend field as 

B ^irh^ < 18 » 

for a given ring size and damping time and, given a max- 
imum wiggler field, the required wiggler length can be 
found from Eqs. i|17|l . At this point, we can calculate 
the other ring parameters. According to the dispersion 
relation in Eq. JBJ, we will need the momentum com- 
paction factor a, the energy spread vq, and the bunch 
rms length a z . Using Eqs. (37), (39) and (40) from Ref. 
PH. the scaling can be written: 



ha — 



6 2irB a 




e 


2?r {Bp) 




R 



cx (1 + F w ) 5/3 ; u cx 

a z cx voy/a . 



B a -\- B W F W 



Examples of the scaling of these parameters versus F w 
can be found in Figs. 11 and 13 of Ref. p^| for an earlier 
version of the NLC damping ring. Using these relations, 
we estimated parameters for NLC ring with F w = and 
F w = 6.5. These parameters are summarized in Table II 
along with the nominal parameters where F w =2.2. 
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FIG. 6: The threshold particle number as a function of the 
CSR wavelength for the NLC main damping ring. The dashed 
curve is for F w — with the vertical dashed straight line 
indicating the corresponding cutoff wavelength. Their cross 
point gives the threshold at the cutoff wavelength. The solid 
curve and vertical solid straight line are for the case of F w = 
2.2, the nominal case in Tabled The dotted curve and vertical 
dotted straight line are for the case of F w = 6.5. The dot- 
dashed curve is for the case of F w — 6.5, but the impedance 
is the dipole CSR impedance alone. 

At this point, we can calculate the impedance. Based 
on our previous studies, we expect the CSR impedance 
from the wiggler to be small compared to the CSR 
impedance of dipoles in the low frequency region of in- 
terest. Hence, the threshold is essentially determined by 
the dipole CSR impedance. According to the impedance 
in Eqs. © and J7J), and the relation of B a and F w in Eq. 
(|15fl . the total impedance should scale as 



Z{k) cx Z D {k)R cx R 1/3 cx (1 + F w ) 1/3 



(20) 



at a fixed wavenumber k. Clearly, this is an increasing 
function with F w however, according to the dispersion 
relation in Eq. ©, the threshold will scale as: 



N t oc 



Z{k) 



cx (I + F™) 13 / 6 , 



(21) 



(19) 



where it was assumed that B a ~ B w , so isq is roughly 
constant. Thus, it is expected that increasing F w will 
increase the threshold significantly. 

In addition, the cutoff wavelength A c cx i?^ 1 / 2 cx 
F^,) -1 / 2 . Hence, as F w increases and B a 



BV 2 cx (1 
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2.2 


6.5 


Dipole field B a /T 


3.8 


1.2 


0.55 


Dipole radius R/m 


1.7 


5.5 


12.0 


Wiggler length L w /m 





46.24 


66 


Momentum compaction a/10 -4 


0.69 


2.95 


14 


Energy rms spread vo/W~ 4 


13 


9.09 


8.4 


Bunch rms length a z /mm 


2.4 


3.6 


5.5 


Cutoff wavelength A c /mm 


8.7 


4.9 


3.3 


Threshold at cutoff At/10 9 


2.0 


5.2 


31.3 



TABLE II: Parameters and results for the NLC main damping ring when F w =0, 2.2 and 6.5 . 



decreases, the cutoff wavelength will decrease. Since, the 
threshold is lowest at long wavelength, the threshold will 
increase as F w is increased due to the shift in the cutoff 
wavelength as well as the change in the ring parameters. 

The above scaling analysis has been compared to nu- 
merical calculations for the NLC damping ring. The 
results are given in Fig. [H] where we plot the thresh- 
old particle number as a function of the CSR instabil- 
ity wavelength for F w = (dashes), 2.2 (solid) and 6.5 
(dots). In addition, the cutoff wavelength for the three 
cases is indicated with vertical lines using the same no- 
tation. Clearly, the threshold increases rapidly as F w 
is increased as expected. The approximately quadratic 
scaling given above is followed as F w is increased from 
the nominal value of 2.2 to 6.5 but the dependence is 
weaker as F w is decreased to because of the change 
in the beam energy spread which was assumed constant. 
Based on these results a new version of the NLC damping 
ring has been designed with an F w 6.5 which has a sig- 
nificantly higher threshold for both this CSR instability 
as well as the conventional microwave instability driven 
by the vacuum chamber impedance |26|. 

VI. DISCUSSION AND CONCLUSION 

In this paper, we studied a possible beam instability 
due to coherent synchrotron radiation from dipole and 
wiggler magnets in a storage ring. This instability may 
be important in many high luminosity factories such as 
DAPHNE, PEP-II, KEK-B, and CESR-C as well as the 
damping rings of future linear colliders. We used the 
results to analyze the stability of the damping ring de- 
signs and found that a possible instability due to the 
CSR will exist in the NLC damping ring design and the 
ATF prototype damping ring at KEK. We then discuss 
optimization of the rings and the pros and cons of in- 
creasing the damping fraction due to the wigglers. This 
led to a redesign for the NLC damping which increases 



the threshold to more than four times the nominal charge 
per bunch. 

As we discussed above, in our calculation, we used 
steady state impedance for both the dipoles and wigglers. 
In the future, it may be important to consider transient 
effects in our calculation. In addition, the shielding cutoff 
of the CSR should be calculated as part of the impedance 
and not treated as an ad hoc limit. This is also important 
because we assumed the same shielding cutoff around the 
ring but in practice it will differ in the bending magnets 
and the wigglers. Finally, the bunch length is the same 
order of the instability wavelength, while in the theory, 
we use a coasting beam model. In the future, the finite 
bunch length should be taken into consideration. All of 
these topics are under study and will be reported in a 
separate paper. 



Acknowledgments 

The authors thank AW. Chao, P. Emma, S.A. Heifets, 
M. Ross, M. Woodley of Stanford Linear Accelerator 
Center, S. Krinsky, J.B. Murphy, J.M. Wang of National 
Synchrotron Light Source, Brookhaven National Labo- 
ratory, and K.J. Kim of Advanced Photon Source, Ar- 
gonne National Laboratory for many discussions. Work 
was supported by the U.S. Department of Energy under 
contract DE-AC03-76SF00515. 



APPENDIX A: GROWTH RATE AT THE FEL 
FREQUENCY AND THE ODD HARMONICS 

In this appendix, we discuss the growth rates calcu- 
lated at the wiggler fundamental frequency, frequently 
refered to as the FEL frequency, and at the odd harmon- 
ics. The wiggler CSR long-range wakefield is 01 
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To the leading order in fc, this gives us the CSR 
impedance as 
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where, is the Heaviside step function, i.e., H(x) = 1 
for x > 0, iJ(0) = 1/2, and = for x < 0. 

At the FEL frequency and the odd harmonics, the 
impedance has a logarithmically divergent imaginary 
part, Im(Z(k)) oc log(|Afc|/fc), with Afc = k— (2n+l)k , 
and the real part is finite. We write Z = R — i L, where 
R > is the resistance and L > is the inductance, 
and L ^S> R. As a first estimate, we neglect the resis- 
tance and write Z = — i L only. In the dispersion rela- 
tion Eq. ©, then w L ^> 1. Now because of the 
factor exp(— p 2 /2), fl is huge compared with p in the p- 
nonvanishing integral region. Hence the integrand could 
be expanded. Keeping the first non-zero term, the inte- 
gral is simplified as 



dp 



p e 



TV2^ 



for large O to get 



a 2 



±- 



LA 



(A3) 



(A4) 



The upper (lower) sign in Eq. (|A4|) refers to the case of 
a positive (negative) rj. For the rings which we study in 
this paper, they all have rj > and we take the upper 
sign. Hence, Im(f2) = 0, i.e., the growth rate will be zero 
and the system is stable at the FEL frequency and all 
the odd harmonics. This agrees with the conventional 
knowledge that a ring is stable with a purely inductive 
impedance and a positive momentum compaction or a 
purely capacitive impedance and a negative momentum 
compaction. 



APPENDIX B: SINGLE WIGGLER 

In this appendix, we study the beam instability in a 
single wiggler. From FEL theory, we would expect in- 
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FIG. 7: The imaginary part of the normalized frequency Q 
as a function of the normalized wavenumber k/ko where ko is 
the on-axis wiggler fundamental radiation wavenumber. The 
solid curve is for the wigglers in the NLC main damping ring 
and the dashed curve is for the TESLA damping ring. 



stability at the wiggler fundamental radiation frequency 
and all the harmonics. In addition, we need to look at a 
microwave-like low-frequency instability. 
In a wiggler, the slippage factor is 



V 



1 + ^ 



(Bl) 



i.e., the slippage factor is negative, therefore, we choose 
the lower (— ) sign in Eq. ©. In Fig. [7] we plot the 
imaginary part of Q as a function of the instability wave 
number for the wigglers in the NLC and TESLA damping 
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rings 0,0]. For the region where Im(f7)< 0, the beam 
is stable. We find that, for the nominal current, the 
system is stable within the entire interesting region where 
the perturbation wavelength A < 3 cm, except at the 
wiggler fundamental radiation wavelength and the third 
harmonic. 

In fact, according to Eq. IjBlfl . we have r\ < and 
we need take the lower sign in Eq. (|A4|) . When we 
approach the fundamental and the third harmonic fre- 
quency, L becomes large, and Eq. (|A4() gives us a large 
Im(f2), i.e., a large growth rate as we find in Fig. [7| How- 
ever, according to Eq. (|A4I) . we will have Im(f2) oc J 1 / 2 , 
where / is the peak current of the electron bunch; while 
in the FEL theory, the growth rate scales as Im(fi) oc J 1 / 3 
[lflj . To correctly calculate the FEL instability from the 
wake impedance approach, the electron density distribu- 
tion function has to be taken at the retarded time in 
the standard Vlasov equation [2(|. This then leads to 
a cubic dispersion relation same as in the FEL theory, 
and results in the cubic root scaling in the growth rate: 



Im(fi) oc I 1 / 3 . 

Nevertheless, here we further use the FEL theory to 
study the instability at the FEL frequency and the third 
harmonic. We use the simulation code TDA [2l| to cal- 
culate the instability growth length at the fundamental 
frequency and a newly developed code (2^ to calculate 
the instability at the third harmonic. The wiggler and 
beam parameters used are listed in Table H] In the NLC 
damping wiggler, the growth length at the FEL frequency 
is about 50 m assuming that the phase errors between 
wiggler sections is small. Fortunately, the separation be- 
tween the 2-meter long wiggler sections will further re- 
duce the gain. In addition, any modulation that is in- 
duced at the fundamental wavelength of ~ 10 /im will be 
smeared out by the energy spread and the momentum 
compaction as the beam is transported around the re- 
mainder of the ring. A similar conclusion is found for 
the TESLA damping ring wigglers and the instability at 
the FEL frequency and the odd harmonics are not ex- 
pected to be a limitation in either case. 
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